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MINOR REVISION 
(October, 1958) 



This edition, Form No. C28-6003, is a minor revision of the preceding 
edition 32-7026-1 but does not obsolete it or 32-7026. The principal 
changes to the former edition are: 



1 1 The statement of the restriction on the value of a subscript expression has been 

revised. 

22 Under Transfer of Control and DOs, the paragraph headed "EXCEPTION" has 

been revised. 

33 The section "READ TAPE" has been revised. 

43 At the beginning of the section "Optimisation of Arithmetic Expressions," sev- 

eral paragraphs have been added on indicating the order of fixed point multi- 
plication and division. 

47 The following sections have been added: "Concept of Relative Constant Defini- 
tion," "Restrictions on Relative Constant Definitions in READ TAPE Lists," 
"Restriction on Relative Constant Definitions in READ DRUM Lists," and 
"Restrictions on Relative Constant Definitions in READ and READ INPUT 
TAPE Lists." 

48 The section "Limits on the Size of the Source Program" has been revised and 
extended. 
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THE FORTRAN SYSTEM 



The IBM Mathematical Formula Translating System Fortran is an automatic 
coding system for the IBM 704 Data Processing System. More precisely, it is a 
704 program which accepts a source program written in the Fortran language, 
closely resembling the ordinary language of mathematics, and which produces 
an object program in 704 machine language, ready to be run on a 704. 

Fortran therefore in effect transforms the 704 into a machine with which 
communication can be made in a language more concise and more familiar than 
the 704 language itself. The result should be a considerable reduction in the 
training required to program, as well as in the time consumed in writing pro- 
grams and eliminating their errors. 
Among the features which characterize the Fortran system are the following. 



Size of Machine 
Required 



The system has been designed to operate on a "small" 704, but to write object 
programs for any 704. (For further details, see the section on Source and Object 
Machines in Chapter 7.) If an object program is produced which is too large 
for the machine on which it is to be run, the programmer must subdivide the 
program. 



Efficiency of the Object programs produced by Fortran will be nearly as efficient as those 
Object Program 



Scope of 
Applicability 



written by good programmers. 

The Fortran language is intended to be capable of expressing any problem of 
numerical computation. In particular, it deals easily with problems containing 



large sets of formulae and many variables, and it permits any variable to have 
up to three independent subscripts. 

However, for problems in which machine words have a logical rather than 
a numerical meaning it is less satisfactory, and it may fail entirely to express 
some such problems. Nevertheless, many logical operations not directly express- 
ible in the Fortran language can be obtained by making use of the provisions 
for incorporating library routines. 



Inclusion of 

Library 

Routines 



Pre-written routines to evaluate any single-valued functions of any number of 
arguments can be made available for incorporation into object programs by 
placing them on the master Fortran tape. 



Provision for Certain statements in the Fortran language cause the object program to be 

Input and Output 

equipped with its necessary input and output programs. Those which deal with 

decimal information include conversion to or from binary, and permit con- 
siderable freedom of format in the external medium. 



Nature of 

Fortran 

Arithmetic 



Arithmetic in the object n ro°ram. will °enerallv be performed with single- 
precision 704 floating point numbers. These numbers provide 27 binary digits 
(about 8 decimal digits) of precision, and may have magnitudes between approxi- 
mately 10~ 38 and 10 38 , and zero. Fixed point arithmetic, but for integers only, 
is also provided. 



CHAPTER 1. GENERAL PROPERTIES 

OF A FQRTRAN so|JRCE pROGRA|V| 



A Fortran source program consists of a sequence of Fortran statements. 
There are 32 different types of statement, which are described in detail in the 
chapters which follow. 



Example of a 
Fortran Program 



The following brief program will serve to illustrate the general appearance and 
some of the properties of a Fortran program. It is shown as coded on a 
standard Fortran coding sheet. 



A FOR 

COMMENT 


6 


FORTRAN STATEMENT 

7 72 


IDENTI- 
FICATION 

73 80 


STATEMENT 

NUMBER 

1 3 


ci 


X 


PROGRAM FOR FINDING THE LARGEST VALUE 




c! 


ATTAINED BY A SET OF NUMBERS 




i 




BIGA = A(lk 








DO 20 I = 2,N 




i 




IF (BIGA - A(D) 10, 20, 20 




i 10 

! 




BIGA = AC I) 




i 20 




CONTINUE 




i 
< i 









This program examines the set of n numbers ai (i=l,. . . ., n) and sets the 
quantity BIGA to the largest value attained in the set. It begins (after a com- 
ment describing the program) by setting BIGA equal to a x . Next the DO 
statement causes the succeeding statements to and including statement 20 to 
be carried out repeatedly, first with i=2, then with i=3, etc., and finally with 
i=n. During each repetition of this loop the IF statement compares BIGA with 
at; if BIGA is less than ai, statement 10, which replaces BIGA by a i5 is executed 
before continuing. 



Keypunching 
the Program 



Each statement is punched on a separate card. If a statement is too long to 
fit on a single card it can be continued over as many as 9 additional continuation 
cards. For each statement the initial card must contain either a zero or a blank 
in column 6; on continuation cards column 6 must not contain a zero or a blank, 
and it should be used to number the continuation cards consecutively from 1 to 9. 
If a statement is too long to fit on a single line of the coding form, the 
programmer can signal to the keypuncher that he has continued on to the next 
line by placing a mark in the column labeled CONTINUATION. 



The order of the statements is governed solely by the order of the cards. 
However, any number less than 2 15 (=32768) may be associated with any 
statement by punching it in columns 1-5 of the initial card bearing that state- 
ment. Thereupon this number becomes the statement number of that statement. 
Statement numbers, besides permitting cross-references within the source pro- 
gram, also help the programmer to correlate the object program with his 
source program. 

Punching the character C in column 1 will cause the card to be ignored by 
Fortran. Such cards may therefore be used to carry comments which will 
appear when the deck is listed. 

Columns 73-80 are not read by Fortran and may be punched with any 
desired identifying information. 

The statements themselves are punched in columns 7-72, both on initial 
and continuation cards. Thus a statement consists of not more than 10 x 66=660 
characters. A table of the admissible characters in Fortran is given in 
Appendix A. 

Blank characters, except in column 6, are simply ignored by Fortran, 
and the programmer may use blanks freely to improve the readability of his 
Fortran listing. 

The actual input to Fortran is either this deck of cards or a BCD tape 
written from it on the 704 peripheral card-to-tape equipment with the standard 
SHARE 80 x 84 board. On the tape an End of File mark after the last card 
is required. 



Preview 

of the Fortran 

Statements 



The 32 types of statement, of which every Fortran program is composed, 
may be classified as follows. 

1. The arithmetic formula, which causes the object program to carry out a 
numerical computation. Chapter 2 discusses the symbols available for refer- 
ring to constants and variables, and Chapter 3 the combining of these into 
arithmetic formulas. 

2. The 15 control statements, which govern the flow of control in the object 
program. These are discussed in Chapter 4. 

3. The 13 input-output statements, which provide the object program with its 
necessary input and output routines. These are discussed in Chapter 5. 

4. Finally, the 3 specification statements, which provide various information 
required or desirable to make the object program efficient, are discussed in 
Chapter 6. 

Chapters 7 and 8, which conclude the manual, give additional detail on various 
topics and examples of Fortran programming. 



CHAPTER 2. THE FORTRAN LANGUAGE: 

Constants, Variables, and Subscripts 



Any programming language must provide for expressing numerical constants and 
variable quantities, Fortran also provides a subscript notation for expressing 
1, 2, or 3-dimensional arrays of variables. 



Constants 



Two types of constant are permissible: fixed point (restricted to integers) and 
floating point (characterized by being written with a decimal point). 

Fixed Point Constants. 



GENERAL FORM 


EXAMPLES 


1 to 5 decimal digits. A preceding + or — 

sign is optional. The magnitude of the constant must 

be less than 32768. 


3 
+1 
-28987 



Any unsigned fixed point constant may be used as a statement number. 
Floating Point Constants. 



GENERAL FORM 


EXAMPLES 


Any number of decimal digits, with 


17. 


a decimal point at the beginning, at the end, 


5,0 


or between two digits. 


-.0003 


A preceding + or — sign is optional. 


5.0E3 (=5.0x103) 


A decimal exponent preceded by an E 


5.0E+3 {= 5.0x103) 


may follow. 


5.0E-7(= 5.0x10-7) 



The magnitude of the number thus expressed must be zero, or must lie between 
the approximate limits of 10 -38 to 10 38 . The number will appear in the object 
program as a normalised single-precision floating point number. 



Variables 



Two types of variable are also permissible: fixed point (restricted to integral 
values) and floating point. Fixed point variables are distinguished by the fact 
that their first character is I, J, K, L, M, or N. ~J" „ ,.^ ^/ 



Fixed Point Variables. 



GENERAL FORM 


EXAMPLES 


1 to 6 alphabetic or numeric characters 
(not special characters) of which the first is 
1, J, K, L, M, or N. 


1 

M2 
JOBNO 



A fixed point variable can assume any integral value whose magnitude is less 
than 32768. However, see the section on Fixed Point Arithmetic in Chapter 7. 

warning. The name of a variable must not be the same as the name of any 
function used in the program after the terminal F of the function name has been 
removed. Also, if a subscripted variable has 4 or more characters in its name, 
the last of these must not be an F. (For the meaning of "function" and "sub- 
scripted" see Chapter 3 and the last section of this chapter.) 

Floating Point Variables. 



GENERAL FORM 


EXAMPLES 


1 to 6 alphabetic or numeric characters 

(not special characters) of which the first is alphabetic 

but not 1, J, K, L F M, or N. 


A 

B7 

DELTA 



A floating point variable can assume any value expressible as a normalised 704 
floating point number; i.e. zero, or with magnitude between approximately 10~ 38 
and 10 38 . 

warning. The restrictions on naming fixed point variables also apply to floating 
point variables. 



Subscripts and 

Subscripted 

Variables 



A variable can be made to represent any member of a 1, 2, or 3-dimensional 
array of quantities by appending to it 1, 2, or 3 subscripts; the variable is then a 
subscripted variable. The subscripts are fixed point quantities whose values 
determine which member of the array is being referred to. 



10 



Subscripts. 



GENERAL FORM 


EXAMPLES 


Let v represent any fixed point variable 


1 


and c (or c') any unsigned fixed point 


3 


constant. Then a subscript is 


MU + 2 


an expression of one of the forms: v 


MU-2 


c 


5*J 


v + cor v— c 


5*J + 2 


c*v 


5 * J— 2 


c*v + c'orc*v— c' 





The symbol * denotes multiplication. The variable v must not itself be sub- 
scripted. 

Subscripted Variables. 



GENERAL F6RM 


EXAMPLES 


A fixed or floating point variable 

followed by parentheses enclosing 1, 2, or 3 

subscripts separated by commas. 


A(l) 
KG) 
BETA(5*J-2, K + 2.L) 



For each variable that appears in subscripted form the size of the array (i.e. the 
maximum values which its subscripts can attain) must be stated in a DIMEN- 
SION statement (see Chapter 6) preceding the first appearance of the variable. 

The value of a subscript expression exclusive of the addend must be greater than 
zero and not greater than the corresponding array dimension. 

Arrangement of Arrays in Storage. 

A 2-dimensional array A will, in the object program, be stored sequentially in 

the order A u , A 2>1 , , A m>1 , A 1>2 , A 2 , 2 , , A m , 2 , , A m , n . Thus 

it is stored "columnwise", with the first of its subscripts varying most rapidly, 
and the last varying least rapidly. The same is true of 3-dimensional arrays, 
i -dimensional arrays are of course simply stored sequentially. All arrays are 
stored backwards in storage; i.e. the above sequence is in the order of decreas- 
ing absolute location. 
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CHAPTER 3. THE FORTRAN LANGUAGE: 

Functions, Expressions, Arithmetic Formulas 
and Function Statements 



Of the various Fortran statements it is the arithmetic formula which defines 
a numerical calculation which the object program is to do. A Fortran arith- 
metic formula resembles very closely a conventional arithmetic formula; it 
consists of the variable to be computed, followed by an = sign, followed by an 
arithmetic expression. For example, the arithmetic formula 

Y = A-SINF(B-C) 
means "replace the value of y by the value of a-sin(b-c)". 



Functions 



As in the above example, a Fortran expression may include the name of a 
function (e.g. the sine function SINF), provided that the routine for evaluating 
the function is either built into Fortran, or is accessible to it as a pre-written 
subroutine in 704 language on the master Fortran tape, or has been defined 
in a function statement (see the last section of this chapter). 



GENERAL FORM 


EXAMPLES 


The name of the function is 4 to 7 alphabetic or numeric 
characters (not special characters), of which the last must 
be F and the first must be alphabetic. Also, the first must 
be X if and only if the value of the function is to be 
fixed point. The name of the function is followed by 
parentheses enclosing the arguments (which may be 
expressions), separated by commas. 


SINF(A+B) 
S0MEF(X,Y) 
SQRTF(SINF(A)) 
XTANF(3.*X) 



Mode of a Function and its Arguments. Consider a function of a single argu- 
ment. It may be desired to state the argument either in fixed or floating point; 
similarly the function itself may be in either of these modes. Thus a function 
of a single argument has 4 possible mode configurations; in general a function 
of n arguments will have 2 n + 1 mode configurations. 

A separate name must be given, and a separate routine must be avail- 
able, for each of the mode configurations which is used. Thus a complete 
set of names for the sine function might be 
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SINOF 


Fixed argument, floating function 


SIN1F 


Floating " , " 


XSINOF 


Fixed " , fixed 


XSIN1F 


Floating " , " 



The X's and F's are compulsory, but the rest of the naming is arbitrary. 

Built-in Functions. The Fortran system has the routines for evaluating certain 
functions built in. The following 18 built-in subroutines are always compiled 
into the object program as open subroutines. 



TYPE OF FUNCTION 


DEFINITION 


NO. OF 
ARGS 


NAME 


MOD 
ARGUMENT 


EOF 

FUNCTION 


Absolute value 


!Arg| 


1 


ABSF 


Floating 


Floating 








XABSF 


Fixed 


Fixed 


Truncation 


Sign of Arg times 


1 


INTF 


Floating 


Floating 




largest integer 




XINTF 


Floating 


Fixed 




^|Arg| 










rvSiTiSiriuSring 


rtigl uituu nigy 


L 


unnr 
iviuur 


Floating 


riuauiig 


(see note below) 






XMODF 


Fixed 


Fixed 


Choosing 


Max (Arg 1? Arg 2 , ..) 


U2 


MAXOF 


Fixed 


Floating 


largest value 






MAX1F 


Floating 


Floating 








XMAXOF 


Fixed 


Fixed 








XMAX1F 


Floating 


Fixed 


Choosing 


Min (Arg 1( Arg 2 , ..) 


£2 


MINOF 


Fixed 


Floating 


smallest value 






MIN1F 


Floating 


Floating 








XMINOF 


Fixed 


Fixed 








XMIN1F 


Floating 


Fixed 


Float 


Float fixed 
number 


1 


FLOATF 


Fixed 


Floating 


Fix 


Same as XINTF 


1 


XFIXF 


Floating 


Fixed 


Transfer 


Sign of Arg 2 


2 


SIGNF 


Floating 


Floating 


of sign 


times Arg, 


2 


XSIGNF 


Fixed 


Fixed 



note. The function MODF (Arg 1? Arg 2 ) is defined as Argx — [Argi/Arg^] 
Arg 2 , where [x] = integral part of x. 



Functions on the Library Tape. Besides the built-in routines, any single- 
valued function of any number of arguments can be made available to the 
programmer by placing the appropriate routine on the master Fortran tape. 
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Any such routine will be compiled into the object program as a closed subrou- 
tine. In the section on Writing Subroutines for the Master Tape in Chapter 7 
are given the specifications which any such routine must meet. 

Expressions An expression is any sequence of constants, variables (subscripted or not sub- 

scripted), and functions, separated by operation symbols, commas, and paren- 
theses so as to form a meaningful mathematical expression. 

However, one special restriction does exist. A Fortran expression may 
be either a fixed or a floating point expression, but it must not be a mixed 
expression. This does not mean that a floating point quantity can not appear 
in a fixed point expression, or vice versa, but rather that a quantity of one 
mode can appear in an expression of the other mode only in certain ways. 
Briefly, a floating point quantity can appear in a fixed point expression only 
as an argument of a function; a fixed point quantity can appear in a floating 
point expression only as an argument of a function, or as a subscript, or as 
an exponent. 

Formal Rules for Forming Expressions. By repeated use of the following 
rules, all permissible expressions may be derived. 

1. Any fixed point (floating point) constant, variable, or subscripted variable 
is an expression of the same mode. Thus 3 and I are fixed point expressions, 
and ALPHA and A(I,J,K) are floating point expressions. 

2. If SOMEF is some function of n variables, and if E, F, , H are a set 

of n expressions of the correct modes for SOMEF, then SOMEF (E, F, 
.... , H) is an expression of the same mode as SOMEF. 

3. If E is an expression, and if its first character is not + or — , then +E and 
— E are expressions of the same mode as E. Thus —A is an expression, but 
-I — A is not. 

4. If E is an expression, then (E) is an expression of the same mode as E. 
Thus (A), ((A)), (((A))), etc. are expressions. 

5. If E and F are expressions of the same mode, and if the first character of 
F is not + or — , then 

E + F 

c c 

1_ — I 

E * F 

E / F 

are expressions of the same mode. Thus A — hB and A/+B are not expres- 
sions. The characters +, — , *, and / denote addition, subtraction, multi- 
plication, and division. 
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6. If E and F are expressions, and F is not floating point unless E is too, and 
the first character of F is not + or — , and neither E nor F is of the form 
A**B, then 

£**p 

is an expression of the same mode as E. Thus A**(B**C) is an expression, 
but I**(B**C) and A**B**C are not. The symbol ** denotes exponentia- 
tion; i.e. A**B means A B . 



Hierarchy of Operations. When the hierarchy of operations in an expression 
is not completely specified by parentheses, then it is understood to be in the 
following order (from innermost operations to outermost) : 

Exponentiation 

Miiltinliratinn and niukinn 

Addition and Subtraction 

For example, the expression 

A + B/C + D**E*F-G 
will be taken to mean 

A + (B/C) + (D E *F)-G. 



Ordering within a Hierarchy. Parentheses which have been omitted from a 
sequence of consecutive multiplications and divisions (or consecutive additions 
and subtractions) will be understood to be grouped from the left. Thus, if . rep- 
resents either * or / (or either + or — ), then 

A.B.C.D.E 
will be taken to mean 

((((A.B).C).D).E) 



Verification of Correct Use of Parentheses. The following procedure is sug- 
gested for checking that the parentheses in a complicated expression correctly 
express the desired operations. 

Label the first open parenthesis "1"; thereafter, working from left to right, 
increase the label by 1 for each open parenthesis and decrease it by 1 for each 
closed parenthesis. The label of the last parenthesis should be 0; the mate of 
an open parenthesis labeled n will be the next parenthesis labeled n-1. 
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Optimisation of Arithmetic Expressions. The efficiency of the object program 
into which an arithmetic expression is translated may depend upon how the 
arithmetic expression is written. The section on Optimisation of Arithmetic 
Expressions in Chapter 7 mentions some of the considerations which affect 
object program efficiency. 



Arithmetic 
Formulas 



GENERAL FORM 


EXAMPLES 


"a=b" where a is a variable (subscripted or 
not subscripted) and b is an expression. 


A(i)=B(l)+SINF(C(l)) 



The = sign in an arithmetic formula has the meaning "is to be replaced by". 
An arithmetic formula is therefore a command to compute the value of the 
right-hand side and to store that value in the storage location designated by the 
left-hand side. 

The result will be stored in fixed or floating point according as the variable 
on the left-hand side is a fixed or floating point variable. 

If the variable on the left is fixed point and the expression on the right is 
floating point, the result will first be computed in floating point and then trun- 
cated and converted to a fixed point integer. Thus, if the result is ±3.569 the 
fixed point number stored will be ±3, not ±4. If the variable on the left is 
floating point and the expression on the right is fixed point, the result will be 
computed in fixed point, truncated to an integer, and converted to floating point. 



Examples of Arithmetic Formulas. 


FORMULA 


MEANING 


A=B 


Store the value of B in A. 


l = B 


Truncate B to an integer, convert to fixed point, 




and store in 1. 


A=l 


Convert 1 to floating point and store in A. 


1 = 1+1 


Add 1 to 1 and store in 1. This example illustrates 




the point that an arithmetic formula is not 




an equation but a command to replace a value. 


A=MAX1F(SINF(B), COSF(B)) 


Replace A by the larger of the quantities 




sinB and cosB. This example illustrates the use 




of a function as an argument of a function. 


A=3.0*B 


Replace A by 3B. 


A=3*B 


Not permitted. The expression is mixed. 


A=I*B 


Not permitted. The expression is mixed. 
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Function 
Statements 



A function may also be denned in the source program itself by means of a func- 
tion statement, and the definition will then persist throughout that one program. 
Thus it is possible to have the convenience of the function notation even for 
functions which are not important enough to deserve a place on the master 
Fortran tape. 





GENERAL FORM EXAMPLES 


"a = b" where a is a function name followed by 
parentheses enclosing its arguments (which must 
be distinct non-subscripted variables) separated 
by commas, and b is an expression (see page 14) 
which does not involve subscripted variables. Any 
functions appearing in b must be built-in, or 
available on the master tape, or already defined 
uy prscsuing lunciion SLaiernents. 


FIRSTF(X) = A*X + B 
SECONDF (X,B) = A*X + B 
THIRDF(D) = FIRSTF(E)/D 
FOURTHF(F,G) = SECONDRF, 
THIRDF(G)) 
FIFTHF(I,A) = 3.0*A**I 
SIXTHF(J) = J + K 
XSiXTHFu) = j + K 



Just as with ordinary functions, the answer will be expressed in fixed or 
floating point according as the name does or does not begin with X. 

The right-hand side of a function statement may be any expression not 
involving subscripted variables that meets the requirements specified for expres- 
sions. In particular, it may involve functions freely, provided that any such 
function, if it is not built-in or available on the master tape, has been defined 
in a preceding function statement. 

As many as desired of the variables appearing in the expression on the 
right-hand side may be stated on the left-hand side to be the arguments of the 
function. Since the arguments are really only dummy variables, their names are 
unimportant (except as indicating fixed or floating point mode) and may even 
be the same as names appearing elsewhere in the program. 

Those variables on the right-hand side which are not stated as arguments 
are treated as parameters. Thus, if FIRSTF is defined in a function statement as 
FIRSTF(X) = A*X + B then a later reference to FIRSTF (Y) will cause 
ay + b, based on the current values of a, b, and y, to be computed. The naming 
of parameters, therefore, must follow the normal rules of uniqueness. 

A function denned by a function statement may be used just as any other 
function. In particular, its arguments may be expressions and may involve 
subscripted variables; thus a reference to FIRSTF (Z + Y(I)) will cause 
a(z + y s ) + b, based on the current values of a, b, y i5 and z, to be computed. 

Functions defined by function statements are always compiled as closed 
subroutines. 

note. All the function statements in a program must precede the first executable 
statement of the program. 
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CHAPTER 4. THE FORTRAN LANGUAGE: 

Control Statements 



Unconditional 
GO TO 



Assigned 
GO TO 



ASSIGN 



The second class of Fortran statements is the set of 15 control statements, 
which enable the programmer to state the flow of his program. 



GENERAL FORM 


EXAMPLES 


"60 TO n" where n is a statement number. 


GO TO 3 



This statement causes transfer of control to the statement with statement 
number n. 



GENERAL FORM 


EXAMPLES 


"GO TO n, foi, rij, , . , , n m )" where n is a 
non-subscripted fixed point variable appearing in a 
previously executed ASSIGN statement, and 
ni, n 2 , , tyn are statement numbers. 


GO TO N, (7, 12, 19) 



This statement causes transfer of control to the statement with statement num- 
ber equal to the value of n last assigned by an ASSIGN statement. The n u n 2 , 
.... , n m are a list of the values which n may be assigned. 
The assigned GO TO is used to obtain a pre-set many-way fork. 

note: When an assigned GO TO exists in the range of a DO, there is a restric- 
tion on the ni, n 2 , . . . , n m . (See the section on DOs in this chapter.) 



GENERAL FORM 


EXAMPLES 


"ASSIGN i TO n" where i is a statement number and 
n is a non-subscripted fixed point variable. 


ASSIGN 12 TO N 



This statement causes a subsequent GO TO n, (n l5 . . , i, . . , n m ) to transfer 
control to the statement whose statement number is i. 

The statement ASSIGN 12 TO N and the arithmetic formula N = 12 are 
not the same. A variable which has been assigned can be used only for an 
assigned GO TO until it is re-established as an ordinary variable. 
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Computed 
GO TO 



IF 



SENSE LIGHT 



IF (SENSE LIGHT) 



GENERAL FORM 


EXAMPLES 


"GO TO (n 1( n 2 , , n m ), i" where 

ni, n 2 , , n m are statement numbers and 

i is a non-subscripted fixed point variable. 


GO TO (30, 40, 50, 60), 1 



If at the time of execution the value of the variable i is j, then control is trans- 
ferred to the statement with statement number nj. Thus, in the example, if I 
has the value 3 at the time of execution, a transfer to statement 50 will occur. 
This statement is used to obtain a computed many-way fork. 



GENERAL FORM 


EXAMPLES 


"IF (a) n,, n 2 , n 3 " where a is any expression 
and n 1( n 2 , n 3 are statement numbers. 


IF (AU.K) -B) 10, 20, 30 



Control is transferred to the statement with statement number n 1? n 2 , or n 3 
according as the value of a is less than, equal to, or greater than zero. 



GENERAL FORM 


1 EXAMPLES 


"SENSE LIGHT i" where i is 0, 1, 2, 3, or 4. 


SENSE LIGHT 3 



If i is 0, all Sense Lights on the 704 console will be turned OFF; otherwise 
Sense Light i will be turned ON. 



GENERAL FORM 


EXAMPLES 


"IF (SENSE LIGHT i) n 1t n 2 " where 

\\i and n 2 are statement numbers and i is 1, 2, 3, or 4. 


IF (SENSE LIGHT 3) 30, 40 



Control is transferred to the statement with statement number n x or n 2 accord- 
ing as Sense Light i is ON or OFF, and the Sense Light is turned OFF. 



IF 

(SENSE SWITCH) 



GENERAL FORM 


EXAMPLES 


"IF (SENSE SWITCH i) n lf n 2 " where n, and n 2 are 
statement numbers and i is 1, 2, 3, 4, 5, or 6. 


IF (SENSE SWITCH 3) 30, 40 
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IF ACCUMULATOR 
OVERFLOW 



IF QUOTIENT 
OVERFLOW 



IF DIVIDE CHECK 



Control is transferred to the statement with statement number n x or n 2 accord- 
ing as Sense Switch i on the 704 console is DOWN or UP. 



GENERAL FORM 


EXAMPLES 


"IF ACCUMULATOR OVERFLOW n,, n 2 " 
where n, and n 2 are statement numbers. 


IF ACCUMULATOR OVERFLOW 30, 40 



Control is transferred to the statement with statement number n x or n 2 accord- 
ing as the Accumulator Overflow trigger of the 704 is ON or OFF, and the 
trigger is turned OFF. 



GENERAL FORM 


EXAMPLES 


"IF QUOTIENT OVERFLOW n,, n 2 " 
where n, and n 2 are statement numbers. 


IF QUOTIENT OVERFLOW 30, 40 



Control is transferred to the statement with statement number n a or n 2 accord- 
ing as the Multiplier-Quotient Overflow trigger of the 704 is ON or OFF, and 
the trigger is turned OFF. 



GENERAL FORM 


EXAMPLES 


"IF DIVIDE CHECK n lf n 2 " where 
n, and n 2 are statement numbers. 


— .,,,■ 

IF DIVIDE CHECK 30, 40 



Control is transferred to the statement with statement number ni or n 2 accord- 
ing as the Divide Check trigger of the 704 is ON or OFF, and the trigger is 
turned OFF. 



PAUSE 



GENERAL FORM 


EXAMPLES 


"PAUSE" or "PAUSE n" where n is an 
unsigned octal fixed point constant. 


PAUSE 
PAUSE 77777 



The machine will HALT, with the octal number n displayed on the 704 con- 
sole in the address field of the storage register. (If n is not stated it is taken 
to be 0.) Pressing the START button causes the program to resume at the 
next Fortran statement. 
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STOP 



GENERAL FORM 


EXAMPLES 


"STOP" or "STOP n" where n is an 
unsigned octal fixed point constant. 


STOP 
STOP 77777 



This statement causes the machine to HALT in such a way that pressing the 
START button has no effect. Therefore, in contrast to the PAUSE, it is used 
where a get-off-the-machine stop, rather than a temporary stop, is desired. The 
octal number n is displayed on the 704 console in the address field of the 
storage register. (If n is not stated it is taken to be 0.) 



DO 



GENERAL FORM 


EXAMPLES 


"DO n i = m ]f m 2 " or "DO n i = m h m 2 , m 3 " 


DO 30 1 = 1, 10 


where n is a statement number, i is a 


DO 30 1 = 1, M, 3 


non-subscripted fixed point variable, and 




m 1( m 2 , m 3 are each either an unsigned fixed point 




constant or a non-subscripted fixed point variable. 




If m 3 is not stated it is taken to be 1. 





The DO statement is a command to "DO the statements which follow, to and 
including the statement with statement number n, repeatedly, the first time with 
i = mi and with i increased by m 3 for each succeeding time; after they have 
been done with i equal to the highest of this sequence of values which does not 
exceed m 2 let control reach the statement following the statement with state- 
ment number n". 

The range of a DO is the set of statements which will be executed re- 
peatedly; it is the sequence of consecutive statements immediately following 
the DO, to and including the statement numbered n. 

The index of a DO is the fixed point variable i, which is controlled by the 
DO in such a way that its value begins at m x and is increased each time by 
m 3 until it is about to exceed m 2 . Throughout the range it is available for com- 
putation, either as an ordinary fixed point variable or as the variable of a 
subscript. During the last execution of the ran^e, the DO is said to be satisfied. 

Suppose, for example, that control has reached statement 10 of the 
program 



10 
11 
12 



DO 11 I = 1, 10 
A(l) = l*N(l) 
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The range of the DO is statement 11, and the index is I. The DO sets I to 1 
and control passes into the range. N(l) is computed, converted to floating 
point, and stored in A(l). Now, since statement 11 is the last statement in the 
range of the DO and the DO is unsatisfied, I is increased to 2 and control 
returns to the beginning of the range, statement 11. 2N(2) is computed and 
stored in A (2) . This continues until statement 1 1 has been executed with 1=10. 
Since the DO is satisfied, control now passes to statement 12. 

DOs within DOs. Among the statements in the range of a DO may be other 
DO statements. When this is so, the following rule must be observed. 

Rule 1. If the range of a DO includes another DO, then all of the statements in the 
range of the latter must also be in the range of the former. 

A set of DOs satisfying this rule is called a nest of DOs. 

Transfer of Control and DOs. Transfers of control by IF-type or GO TO- 
type statements are subject to the following rule. 

Rule 2. No transfer is permitted into the range of any DO from outside its range. 
Thus, in the configuration below, 1, 2, and 3 are permitted transfers, but 
4, 5, and 6 are not. 




exception. There is one situation in which control can be transferred into the 
range of a DO from outside its range. Suppose control is in the range of the 
innermost DO of a nest of DOs which are completely nested: i.e., for every 
pair of DOs in the nest, one contains the other. Suppose also that control is 
transferred to a section of program, completely outside the nest to which those 
DOs belong, which makes no change in any of the indexes or indexing parame- 
ters (m's) in the nest. Then after the execution of this section of program, 
control can be transferred back to the range of the same innermost DO from 
which it originally came. This provision makes it possible to exit temporarily 
from the range of some DOs to execute a subroutine. An example is given in 
Chapter 8. 
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Restriction on Assigned GO TO's in the Range of a DO. When an assigned 
GO TO is in the range of a DO, the statements in the nest to which it may trans- 
fer must all be in the exclusive range of a single DO, i.e., among those state- 
ments in the range of a DO which are not in the range of any DO in its range. 

Preservation of Index Values. When control leaves the range of a DO in the 
ordinary way (i.e. by the DO becoming satisfied and control passing on to the 
next statement after the range) the exit is said to be a normal exit. After a 
normal exit from a DO occurs, the value of the index controlled by that DO 
is not defined, and the index can not be used again until it is redefined. (See, 
however, the section on Further Details about DO Statements in Chapter 7.) 
However, if exit occurs by a transfer out of the range, the current value 
of the index remains available for any subsequent use. If exit occurs by a trans- 

ici wiiivij. is m uic lcuigoa ui acvciai t/vjs, uic wuiicin vaiucs ui an luc liiucAca 

controlled by those DOs are preserved for any subsequent use. 

Restriction on Calculations in the Range of a DO. Almost every type of 
calculation is permitted in the range of a DO. Only one type of statement is 

"Vt ^VlilllUVU) XXC4.XXXWXJT W-AXJT OUllVlUVlll VfillVil XWU.WJ_LJ.XWt3 IXXW VU1UV WX IXXW 111UUA 

or of any of the indexing parameters (m's). In other words, the indexing of a 
DO loop must be completely set before the range is entered. 

The first statement in the range of a DO must not be non-executable. 



CONTINUE 



GENERAL FORM 


EXAMPLES 


"CONTINUE" 


C0NTINUE 



CONTINUE is a dummy statement which gives rise to no instructions in the 
object program. Its most frequent use is as the last statement in the range of a 
DO, where it provides a statement number which can be referred to in transfers 
which are desired to become, in the object program, transfers to the indexing 
instuctions at the end of the loop. 

As an example of a program which requires a CONTINUE, consider the 



iaoie SCarun program 



10 DO 12 I = 1, 100 

11 IF{ARG— VALUEd) ) 12, 20, 12 

12 CONTINUE 
13 
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This program will examine the 100-entry VALUE table until it finds an entry 
which equals ARG, whereupon it will exit to statement 20 with the successful 
value of I available for fixed point use; if no entry in the table equals ARG 
a normal exit to statement 13 will occur. The program 

10 DO 11 1 = 1,100 

11 IF(ARG-VALUE(1) ) 11, 20, 11 

12 

would not work since, as stated in the next section, DO-sequencing does not 
occur if the last statement in the range of a DO is a transfer. 



Summary of 

Fortran 

Sequencing 



The precise laws which govern the order in which the statements of a Fortran 
program will be executed, and which have been left unstated till now, may be 
stated as follows. 

1. Control begins at the first executable statement. 

2. If control is in a statement S, then control will next go to the statement dic- 
tated by the normal sequencing properties of S. (The normal sequencing 
for each type of Fortran statement is given in Appendix B.) 

3. EXCEPTION. If, however, S is the last statement in the range of one or 
more DOs which are not yet satisfied, and if S is not a transfer (IF-type or 
GO TO-type statement), then the normal sequencing of S is ignored and 
DO-sequencing occurs; i.e. control will next go to the first statement of the 
range of the nearest of the unsatisfied DOs, and the index of that DO will 
be raised. 

note. The statements FORMAT, DIMENSION, EQUIVALENCE, and FRE- 
QUENCY, which are discussed in the next two chapters, are non-executable 
statements, and in any question of sequencing are simply to be ignored. 

If the last executable statement in the source program is not a STOP or 
an IF-type or GO TO-type transfer, then the object program will be compiled 
as if the pseudo-statement PRESS LOAD CARDS button had followed the last 
executable statement. 



warning. Every executable statement in the source program (except the first) 
must have some path of control leading to it; otherwise errors will occur in the 
compilation of the object program. 
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CHAPTER 5. THE FORTRAN LANGUAGE: 



There are 13 Fortran statements available for specifying the transmission of 
information, during the running of the object program, between core storage 
on the one hand and tapes, drums, and attached card reader, card punch, and 
printer on the other hand. These 13 statements may be grouped as follows. 

1. Five statements (READ, PUNCH, PRINT, READ INPUT TAPE, and 
WRITE OUTPUT TAPE) call for the transmission of a list of quantities 
between cores and an external storage medium — cards, printed sheet, or 
tape — in which information is expressed in Hollerith punch, decimal print, 
or BCD code, respectively. 

2. One statement (FORMAT) is a non-executed statement which defines the 
information format in the external medium for the above 5 statements. 

3. Four more statements (READ TAPE, READ DRUM, WRITE TAPE, and 
WRITE DRUM) call for ordinary binary transmission of a list of quantities. 

4. Finally, there are 3 statements (END FILE, REWIND, and BACKSPACE) 
for manipulating tapes. 

This chapter will first discuss the formation of a list of quantities for trans- 
mission, since such a list must appear in each of the 9 statements in groups 
1 and 3 above. Next the method of writing a FORMAT statement will be 
described, and the format which input data to the object program must have. 
Finally, the statements in groups 1, 3, and 4 will be discussed. 



Specifying Each of the 9 statements which call for the transmission of information includes 

a list of the quantities to be transmitted. This list is ordered, and its order must 



Quantities 



be the same as the order in which the words of information exist (for input), 
or will exist (for output), in the external medium. 

The formation and meaning of a list is best described by an example. 
A, B(3), (C(l), D(I,K), I = 1,10), ((E(I,J), I = 1,10,1), F(J,3), J = 1,K) 

Suppose that this list is used with an output statement. Then the information 
will be written in the external medium in the order 
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A, B0), C(l), D(1,K), C(2), D(2,K), , C(10), D(10,K), 

E(l,l), E(2,l), , E(10,l), F(l,3), 

E(1,2),E(2,2), ,E(10,2),F(2,3), ,F(K,3). 

Similarly, if this list were used with an input statement, the successive words, 
as they were read from the external medium, would be placed into the sequence 
of storage locations just given. 

Thus the list reads from left to right and with repetition of variables en- 
closed within parentheses. Only variables, and not constants, may be listed. 
The repetition is exactly that of DO-repetition, as if each open parenthesis 
(except subscripting parentheses) were a DO, with indexing given immediately 
before the mated closed parenthesis, and with range extending -up to that index- 
ing information. The order of the above list is the same as of the "program" 



1 


A 


2 


B(3) 


3 


0051 = 1,10 


4 


CO) 


5 


D(I,K) 


6 


DO 9 J = 1, K 


7 


D08l = l,10,l 


8 


E(I,J) 


9 


F(J,3) 



Notice that indexing information, as in DOs, consists of 3 constants or 
fixed point variables, and that the last of these may be omitted, in which case 
it is taken to be 1. 

For a list of the form K, A(K) or K, (A(I), I = 1,K), where an index 
or indexing parameter itself appears earlier in the list of an input statement, 
the indexing will be carried out with the newly read-in value. 

When it is desired to transmit an entire array, and in the inverse of its 
natural order, i.e., in the order A,„. n , . . . , A mt ->, . . . , A 2 .j, A 1>2 , A„ U) . . . , 
A 2 ,i, Ai,i, etc. (see the section on Arrangement of Arrays in Storage in Chapter 
2), then an abbreviated notation may be used; only the name of the array need 
be given and the indexing information may be omitted. Fortran will look to see 
if a DIMENSION statement (see Chapter 6) has been made about that name; 
if it has, indexing information to transmit the entire array in the inverse of its 
natural order will be supplied automatically, while if it has not, only a single 
variable will be transmitted. Thus, in the example, the entire A-array will be 
transmitted, including the case where the array consists of a single quantity. 
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FORMAT 



warning. The information given in this section applies in its full generality 
only to lists which are given with the 5 decimal statements of group 1 , For the 
binary statements WRITE DRUM and READ DRUM only the abbreviated 
notation mentioned immediately above is permitted; the restrictions which affect 
lists for WRITE TAPE and READ TAPE are discussed in the section on Lists 
for Binary Tape Operations in Chapter 7. 



GENERAL FORM 


EXAMPLES 


"FORMAT (Specification)" where Specification 
is as described below. 


FORMAT (I2/CE12.4, F10.4)) 



The 5 decimal input-output statements of group 1 contain, in addition to the 
list of quantities to be transmitted, the statement number of a FORMAT state- 
ment, which describes the information format which exists, or is to be produced, 
in the external medium. It also specifies the type of conversion between binary 
and decimal which is to be performed. FORMAT statements are not executed, 
their function being merely to supply information in the object program. There- 

The Specification. For the sake of definiteness the details of writing a FORMAT 
Specification are given for use with a PRINT statement. However, the descrip- 
tion is valid for any case simply by generalizing the concept of "printed line" 
to that of unit record in the external medium. Thus a unit record may be 

1. A printed line with a maximum of 120 characters. 

2. A punched card with a maximum of 72 characters. 

3. A BCD tape record with a maximum of 120 characters. 

Three basic types of decimal-to-binary or binary-to-decimal conversion are 
available: 



TO OR FROM 



EXTERNAL 



Floating point variable 
Floating point variable 
Fixed point variable 



Floating point decimal 
Fixed point decimal 
Decimal integer 



The FORMAT specification describes the line to be printed by giving, for each 
field in the line (from left to right, beginning with the first type wheel) : 
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1. The type of conversion (E, F, or I) to be used; 

2. The width (w) of the field; and 

3. For E- and F-conversion, the number of places (d) after the decimal point 
that are to be rounded and printed. If d is not less than 10 it is treated 
mod 10. 

These basic field specifications are given in the forms 

Iw, Ew.d, and Fw.d 

with the specification for successive fields separated by commas. Thus the state- 
ment FORMAT (12, E12.4, F10.4) might give the line 

27 -0.9321E 02 -0.0076 

As in this example the field widths may be made greater than necessary, so as 
to ensure spacing blanks between numbers. In this case there is 1 blank follow- 
ing the 27, 1 blank automatically supplied after the E, and 3 blanks after the 
02. Within each field the printed output is always pushed to the extreme right. 
It may be desired to print n successive fields in the same fashion. This may 
be done by giving n before the E, F, or I. Thus the statement FORMAT (12, 
3E12.4) might give 

27 -0.9321E 02 -0J580E-02 0.5536E 00 

To permit the repetition of groups of field specifications a limited parenthetical 
expression is permitted. Thus FORMAT (2(F10.6, E10.2), 14) is equivalent 
to FORMAT (F10.6, E10.2, F10.6, E10.2, 14). No provision is made for 
parentheses within parentheses. 

To permit more general use of F-conversion, a scale factor followed by the 
letter P may precede the specification. The scale factor is so defined that 

Printed number = Internal number X io Scalefactor . 

Thus the statement FORMAT (12, 1P3F1.1.3), used with the data of the 
preceding example, would give 

27 -932.096 -0.076 5.536 

while FORMAT (12, -1P3F11.3) would give 

27 -9.321 -0.001 0.055 

A positive scale factor may also be used with E-conversion to increase the num- 
ber and decrease the exponent. Thus FORMAT (12, 1P3E12.4) would give 
with the same data 

27 -9.3210E 01 -7.5804E-03 5.5361E-01 

28 



note. The scale factor is assumed to be zero if no other value has been given. 
However, once a value has been given it will hold for all E- and F-conversions 
following the scale factor within the same FORMAT statement. This applies to 
both single-record and multi-record formats (see below). Once a scale factor 
has been given, a subsequent scale factor of zero in the same FORMAT state- 
ment must be specified by OP. Scale factors have no effect on I-conversion. 

Hollerith Fields. A field may be designated as a Hollerith field, in which case 
English text will be printed in it. The field width, followed by the desired 
characters, should appear in the appropriate place in the specification. For 
example, FORMAT (3HXY= F8.3, 4H Z = F6.2, 7H W/AF = F7.3) 
would give with the same data 

XY= -93.210 1= -0.01 W/AF= 0.554 
Notice that any Hollerith characters, including blanks, may be printed. This is 
the sole exception to the statement in Chapter 1 that Fortran ignores blanks. 

It is possible to print Hollerith information only, by giving no list with the 
input-output statement and setting up no I, E, or F fields in the FORMAT 
statement. 

Consider a Hollerith field in a FORMAT statement at the time of execu- 
tion of the object program. If the FORMAT statement is being used with an 
input statement, the Hollerith text in the FORMAT statement will be replaced 
with whatever text is read from the corresponding field in the external medium. 
When the FORMAT statement is used with an output statement, whatever 
text is currently in the FORMAT statement will be outputted. Thus text can 
be originated either at source time or at object time by not using, or using, the 
FORMAT statement with an input statement. 

Multi-Record Formats. To deal with a block of printing a FORMAT speci- 
fication may have several different line formats, separated by the slash /.. Thus 
FORMAT (3F9.2, 2F10.4/8E14.5) would specify a block in which lines 

1, 3, 5, have format 3F9.2, 2F10.4 and lines 2, 4, 6, have format 

8E14.5. 

If a block format is desired havin ff the first two lines of some special 
formats and all the remaining lines of another format, its last line should be en- 
closed in parentheses; e.g. FORMAT (12, 3E12.4/2F10.3, 3F9.4/( 10F12.4) ). 

In general, if there are items in the list still to be transmitted after the 
format specification has been completely used, the format repeats from the last 
open parenthesis or (if no parentheses are present) from the beginning. 
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Blank lines may be introduced into a block by omitting format information; 
thus // and /// will give 1 and 2 blank lines respectively. 

Ending a FORMAT Statement. A decimal input-output operation will be 
brought to an end when and only when a numerical field or the end of the 
FORMAT statement is encountered and there are no items remaining in the 
list of the input-output statement itself. 

Carriage Control. The peripheral printer can operate in 3 modes: single space, 
double space, and Program Control. Under Program Control, which gives the 
greatest flexibility, the first character of each BCD record controls the carriage, 
and that character is not printed. The control characters and their effects are 

Blank Single space before printing 

Double space before printing 

+ No space before printing 

1-9 Skip to channels 1-9 

J-R Short skip to channels 1-9 

Thus a FORMAT Specification for WRITE OUTPUT TAPE for printing with 
Program Control will usually begin with IH followed by the appropriate control 
character. The same is true of PRINT since in Fortran printing on the attached 
printer simulates Program Control printing on the peripheral printer. 

Data Input Decimal data to be read by a READ or READ INPUT TAPE at the time of 

tothe execution of the object program must be in essentially the same format as 

Object Program 

given in the examples of the preceding section. Thus a card to be read with 

FORMAT (12, El 2.4, F10.4) might be punched 

27 -0.9321E 02 -0.0076 
Within each field all information must be pushed to the extreme right. Positive 
signs may be indicated either by a blank or a + ; - signs may be punched with 
an 11 -punch or an 8-4 punch (see Appendix A). Blanks in numeric fields are 
regarded as zeroes. Numbers for E- and F-conversion may contain any number 
of digits, but only 8 digits of accuracy will be retained. Numbers for I-conversion 
will be treated mod 32768. To permit economy in punching certain relaxations 
in input data format are permitted. 

1. Numbers for E-conversion need not have 4 columns devoted to the ex- 
ponent field. The start of the exponent field must be marked by an E, or 
if that is omitted, by a + or - (not a blank). Thus E2, E02, +2, + 02, 
E 02, and E+02 are all permissible exponent fields. 
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READ 



2. Numbers for E- or F-conversion need not have the decimal point punched. 
If it is not punched the FORMAT Specification sets its effective position; 
for example, —09321+2 with E12.4 will be treated as if the decimal point 
had been punched 4 places before the start of the exponent field, that is 
between the and the 9. If the decimal point is punched, its position over- 
rides the value of d given in the FORMAT Specification. 



GENERAL FORM 


EXAMPLES 


"READ n, List" where n is the 

statement number of a FORMAT statement, 

and List is as previously described. 


READ 30, K, A(K) 



The READ statement causes the object program to read cards from the attached 
card reader. Record after record (card after card) is read until the complete 
list has been brought in, converted, and stored in the locations specified by the 
list. The FORMAT statement describes the arrangement of information on the 
cards and the type of conversion to be done. 

If an End of File is encountered (the program attempts to read a card 
and finds that there is none in the card reader) the object program HALTS. 
Placing more cards in the card reader and pressing the START button causes 
the program to continue the reading from the point in the list which it had 
reached. 

A partial check is made in the object program for incorrectly punched 
columns. Such a column causes a HALT. Pressing the START button causes 
the faulty column to be treated as a zero, and the program to continue. 



READ INPUT 
TAPE 



GENERAL FORM 


EXAMPLES 


"READ INPUT TAPE i, n, List" where i is 


READ INPUT TAPE 3, 30, K, A(K) 


an unsigned fixed point constant 


READ INPUT TAPE 1, 30, K,A(K) 


between 1 and 10 inclusive or a fixed point 




variable, n is the statement number 




of a FORMAT statement, and List is as 




picviuuoijr ucooi iucu. 





The READ INPUT TAPE statement causes the object program to read BCD 
information from tape unit i, where i = 1, 2, . . ., 10. Record after record 
is brought in, in accordance with the FORMAT statement, until the complete 
list has been placed in storage. 
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PUNCH 



PRINT 



An End of File causes a HALT in the object program. Pressing the START 
button causes the program to continue the reading from the point in the list 
which it had reached. 

The object program redundancy checks the tape reading. If a record fails 
twice the program HALTS. Pressing the START button causes the information 
read on the second attempt to be accepted and the program to continue. 



GENERAL FORM 


EXAMPLES 


"PUNCH n, List" where n is the statement number of a 
FORMAT statement and List is as previously described. 


PUNCH 30, (A(J), J = 1, 10) 



The PUNCH statement causes the object program to punch cards on the at- 
tached card punch. Card after card is punched in accordance with the FORMAT 
statement until the complete list has been punched. 

No checking is done, and there are no HALTS in the object program. 



GENERAL FORM 


EXAMPLES 


"PRINT n, List" where n is the statement number of a 
FORMAT statement and List is as previously described. 


PRINT 30, (A(J), J = 1, 10) 



The PRINT statement causes the object program to print on the attached 
printer. Line after line is printed in accordance with the FORMAT statement 
until the complete list has been printed. 

The printing is echo checked. A printing error so detected causes the 
object program to HALT. Pressing the START button causes the program to 
continue. Pressing the RESET button and then the START button causes the 
line to be printed again and the program to continue. 



WRITE 
OUTPUT TAPE 



GENERAL FORM 


EXAMPLES 


"WRITE OUTPUT TAPE i, n, List" 
where i is an unsigned fixed 
point constant between 1 and 10 inclusive 
or a fixed point variable, n is the 
statement number of a FORMAT statement, 
and List is as previously described. 


WRITE OUTPUT TAPE 3, 30, (AU), J = 1, 10) 
WRITE OUTPUT TAPE 1, 30, (AU), J = 1, 10) 
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READ TAPE 



The WRITE OUTPUT TAPE statement causes the object program to write 
BCD information on tape unit i, where i = 1, 2, . . ., 10. Record after record 
is written in accordance with the FORMAT statement until the complete list 
has been written. No End of File is written after the last record. 

No checking is done, and there are no HALTS in the object program. 



GENERAL FORM 


EXAMPLES 


"READ TAPE i, List" where i is an 
unsigned fixed point constant between 
1 and 10 inclusive or a fixed 
point variable, and List is as described 
in Chapter 7. 


READ TAPE 3, (AU), J = 1, 10) 
READ TAPE 1, (AU), J = 1, 10) 



The READ TAPE statement causes the object program to read binary informa- 
tion from tape unit i, where i = 1, 2, . . . ., 10. Only one record is read, and 
it will be completely read only if the list contains as many words as the record. 
The tape, however, always moves all the way to the next record. 

The first word of every binary record is copied into an erasable location. 
This copy serves as an End-of-File test. This means that the binary record 
must always contain at least one more word than the binary list. If the record 
has been written by a Fortran compiled program, the extra first word is 
automatically produced. If it is not, care should be taken to write it into the 
record. 

The READ TAPE list should not be longer than the record. If it is, the 
program may stop with a Read-Write check or it may go on to the next 
Fortran statement. This depends on the size of the list remaining and the 
extent of indexing instructions required by the list. 

If an End-of-File or empty record is encountered while doing Fortran 
reading, the test that occurs with the first word takes care of it. In the End-of- 
File case, this test causes a halt in the object program. Pressing the START 
button causes a repetition of this test on the next file. In the Empty Record 
case, the repetition of the test on the next record occurs automatically. (An 
empty record may result from a hand-coded program; for example, where two 
successive write-select tape instructions occur.) 

When, by this process, a real record of information is encountered, it is 
read in according to the list of the READ TAPE instruction. If the instruction 
has no list, the record is spaced over and the next Fortran instruction is 
executed. 
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READ DRUM 



For reasons of timing there are limitations on the complexity of the list. See 
the section on Lists for Binary Tape Operations in Chapter 7. 

The object program redundancy checks the tape reading. (The longitudinal 
check is applied only if the whole record is read.) If a record fails twice the 
program HALTS. Pressing the START button causes the information read on 
the second attempt to be accepted and the program to continue. 



GENERAL FORM 


EXAMPLES 


"READ DRUM i, j, List" where 
i and j are each either an unsigned 
fixed point constant or a fixed 
point variabie, with the vaiue of i 
between 1 and 8 inclusive, and 
List is as described below. 


READ DRUM 2, 1000, A, B, C, D 
READ DRUM 1, J, A, B, C, D 



The READ DRUM statement causes the object program to read words of 
binary information from consecutive locations on drum i, where i = 1, 2, . . . . , 
8, beginning with the word in drum location j, where j = 0, 1, . . . ., 2047. 
(If j ^2048 it is interpreted mod 2048.) Reading continues until the complete 
list has been read in. 

For reasons of timing there are stringent limitations on the complexity of 
the list. In fact, the list can employ only the abbreviated notation described 
earlier in this chapter; it may consist only of variables without subscripts, as 
A, B, C, D, . . . Those variables which are simple will be read into storage in 
the ordinary way; those which are arrays will be read with indexing obtained 
from their DIMENSION statements (see Chapter 6). Thus with READ DRUM 
the full array must be read in, and in the inverse of natural order. 



WRITE TAPE 



GENERAL FORM 


EXAMPLES 


"WRITE TAPE i, List" where i is an 
unsigned fixed point constant between 1 and 10 
inclusive or a fixed point variable, 
and List is as described in Chapter 7. 


WRITE TAPE 3, (AU), J = 1, 10) 
WRITE TAPE 1, (AU), J = 1, 10) 



The WRITE TAPE statement causes the object program to write binary infor- 
mation on tape unit i, where i = 1, 2, . . ., 10. Only one record is written; its 
length will be that of the list, except that an extra first word is written to serve 
in the End-of-File test during reading. 
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WRITE DRUM 



END FILE 



REWIND 



For reasons of timing, there are limitations on the complexity of the list. See 
the section on Lists for Binary Tape Operations in Chapter 7. 

No checking is done, and there are no HALTS in the object program. 



GENERAL FORM 


EXAMPLES 


"WRITE DRUM i, j, List" where' i and j 
are each either an unsigned fixed 
point constant or a fixed point variable, 
with the value of i between 1 and 8 inclusive, 
and List is as described for READ DRUM. 


WRITE DRUM 2, 1000, A, B, C, D 
WRITE DRUM 1, J, A, B, C, D 



The WRITE DRUM statement causes the object program to write words of 
binary information into consecutive locations on drum i, where i = 1, 2, . . . ., 
8, beginning with drum location j,' where j = 0, 1, . . . ., 2047. (If j ^2048 
it is interpreted mod 2048.) Writing continues until the complete list has been 
written. 

The list is subject to the same restrictions as for READ DRUM. 

No checking is done and there are no HALTS in the object program. 



GENERAL FORM 


EXAMPLES 


"END FILE i" where i is an unsigned fixed 
point constant between 1 and 10 inclusive 
or a fixed point variable. 


END FILE 3 
END FILE I 



The END FILE statement causes the object program to write End of File 
on tape unit i, where i = 1, 2, . . . . , 10. 



GENERAL FORM 


EXAMPLES 


"REWIND i" where i is an unsigned fixed 
point constant between 1 and 10 inclusive 
or a fixed point variable. 


REWIND 3 
REWIND 1 



The REWIND statement causes the object program to rewind tape unit i, where 
i = 1,2, ,10. 
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BACKSPACE 



GENERAL FORM 


EXAMPLES 


"BACKSPACE i" where i is an unsigned fixed 
point constant between 1 and 10 inclusive 
or a fixed point variable. 


BACKSPACE 3 
BACKSPACE 1 



The BACKSPACE statement causes the object program to backspace tape 
unit i by one record, where i = 1, 2, . . . . , 10. 

Error Halts The several HALTS which can occur during input or output operation in the 

object program can be identified by the contents of the storage register on the 
704 console. 
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CHAPTER 6. THE FORTRAN LANGUAGE: 

Specification Statements 



The last class of Fortran statements is the set of 3 specification statements 
DIMENSION, EQUIVALENCE, and FREQUENCY. These are statements 
which are not executed, but which furnish information for use by Fortran to 
make the object program efficient. 



DIMENSION 



GENERAL FORM 


EXAMPLES 


"DIMENSION v, v, v, ..." where each 
v is a variable subscripted with 1, 2, or 3 
unsigned fixed point constants. 
Any number of v's may be given. 


DIMENSION A(10), B(5, 15), CO, 4, 5) 



The DIMENSION statement provides the information necessary to allocate 
storage in the object program for arrays of quantities. 

Every variable which appears in the program in subscripted form must 
appear in a DIMENSION statement, and the DIMENSION statement must 
precede the first appearance of the variable. In the DIMENSION statement are 
given the desired dimensions of the array; in the executed program the sub- 
scripts of that variable must never take on values larger than those dimensions. 

Thus the example states that B is a 2-dimensional array and that the sub- 
scripts of B will never exceed 5 and 15; it causes 75 words of storage to be 
set aside for the B array. 

A single DIMENSION statement may be used to dimension any number 
of arrays. 
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EQUIVALENCE 



GENERAL FORM 


EXAMPLES 


"EQUIVALENCE (a, b, c, ...), (d, e, f, ...), .." 
where a, b, c, d, e, f, ... are variables 
optionally followed by a single unsigned 
fixed point constant in parentheses. 


EQUIVALENCE (A, B(l), C(5) ), (D(17), EG) ) 



The EQUIVALENCE statement enables the programmer, if he wishes, to con- 
trol the allocation of data storage in the object program. In particular, it permits 
him to economise on data storage requirements by causing storage locations 
to be shared by two or more quantities, when the logic of his program permits. 

An EQUIVALENCE statement may be placed anywhere in the source 
program. Each pair of parentheses encloses the names of two or more quanti- 
ties whose storage locations are to be made the same in the object program; 
any number of equivalences (pairs of parentheses) may be given. 

In an EQUIVALENCE statement the meaning of C(5), for example, 
is "the 4th storage location in the object program after the cell containing 
C, or (if C is an array) after C(l) or C(l,l) or C(l,l,l)". In general A(p) 
is defined for p ^ 1 and means the p-l th location after A or the beginning of the 
A array; i.e. the p th location in the array. If p is not given, it is taken to be 1 . 

Thus the example statement causes A, B, and C (or the beginnings of 
the A, B, and C arrays) to be so placed in storage that the location containing 
A, the location containing B, and the 4th location after that containing C, are 
the same location. Similarly, it causes the 16th location after D and the 2nd 
after E both to be another location. 

A quantity or array which does not appear in any EQUIVALENCE state- 
ment will have storage exclusively to itself. 

Locations can be shared only among variables, not among constants. 

The sharing of storage locations cannot be planned safely without a 
knowledge of which Fortran statements, when executed in the object pro- 
gram, will cause a new value to be stored in a storage location. There are 7 
such statements. 

1. Execution of an arithmetic formula stores a new value of the variable on its 
left-hand side. 

2. Execution of an ASSIGN i TO n stores a new value in n. 
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3. Execution of a DO will in general store a new value of the index. (It will 
not always do so, however; see the section on Further Details about DO 
Statements in Chapter 7.) 

4. Execution of a READ, READ INPUT TAPE, READ TAPE, or READ 
DRUM stores new values of the variables listed. 



FREQUENCY 



GENERAL FORM 


EXAMPLES 


"FREQUENCY n(i, j, ...), m(k, 1, ...),..." 
where n, m, ... are statement numbers and 
i, j, k, 1, . . . are unsigned fixed point constants. 


FREQUENCY 30(1, 2, 1), 
40(11), 50(1, 7, 1, 1) 



The FREQUENCY statement permits the programmer to give his estimate, for 
each branch-point of control, of the frequencies with which the several branches 
will actually be executed in the object program. This information is used to 
optimise the use of index registers in the object program. 

A FREQUENCY statement may be placed anywhere in the object pro- 
gram, and may be used to give the frequency information about any number 
of branch-points. For each branch-point the information consists of the state- 
ment number of the statement causing the branch, followed by parentheses en- 
closing the estimated frequencies separated by commas. 

Consider the example. This might be a FREQUENCY statement in a 
program in which statement 30 is an IF, 40 is a DO, and 50 is a computed 
GO TO. The programmer estimates that the argument of the IF is as likely to 
be zero as non-zero, and when it is non-zero it is as likely to be negative as 
positive. The DO statement at 40 is presumably one for which at least one of 
the indexing parameters (m's) is not a constant but a variable, so that the num- 
ber of times the loop must be executed to make a normal exit is not known in 
advance; the programmer here estimates that 11 is a good average for that 
number. The computed GO TO at 50 is estimated to transfer to its four branches 
with frequencies 1, 7, 1, 1. 

All frequency estimates, except those about DOs, are relative; thus they 
can be multiplied by any constant. The example statement, for instance, could 
equally well be given as FREQUENCY 30(2,4,2), 40(11), 50(3,21,3,3). A 
frequency may be estimated as 0; this will be taken to mean that the frequency 
is very small. 

The following table lists the 7 types of statement about which frequency 
information may be given. 
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TYPE 


NO OF 
FREQS 


REMARKS 


Computed GO TO 

IF 

IF (SENSE SWITCH) 

IF ACCUMULATOR OVERFLOW 

IF QUOTIENT OVERFLOW 

IF DIVIDE CHECK 

DO 


^2 
3 
2 
2 
2 
2 
1 


Order of frequencies: same as order of branches 
To be given only when m,, m 2 , or m 3 is variable. 



It is not necessary to give frequency information about any branch-point. If 
none is given, it will be taken that the probabilities of all branches are equal. 

A frequency estimate concerning a DO will be ignored except when at 
least one of the indexing parameters of that DO is variable. Moreover, the fre- 
quency estimate should be based only on the expected values of those para- 
meters; in other words, even if the range of the DO contains IFs or GO TOs 
which may transfer outside the range, the frequency estimate should be the 
number of times the range must be executed to cause a normal exit. 

A DO for which the indexing parameters are variable and for which no 
FREQUENCY statement is given will be treated as if a frequency of 5 had 
been estimated. 
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CHAPTER 7. MISCELLANEOUS DETAILS ABOUT FORTRAN 



Source and 
Object Machines 



The source machine is the 704 on which the source program is translated into the 
object program. The object machine is that on which the object program is run. 

The source machine must be at least as large as a "small" 704; i.e. a 704 pos- 
sessing 4096 words of core storage, floating point instructions, CP A (copy and add 
logical) instruction, 1 drum unit, 4 tape units, attached card punch, attached or 
peripheral card reader, and attached or peripheral printer. 

The object machine may be of any size. Among the information produced by 
Fortran is a count of the storage locations required by the object program, from 
which it can be easily decided whether the object program is too large for any given 
object machine. 



Arrangement 

of the 

Object Program 



The instructions and constants of the object program begin in lower memory and 
extend upwards. Data and other storage locations required for the operation of the 
program begin at location 11111% and extend downwards. Thus these latter locations 
are always at the top of memory, regardless of the size of the object machine. 

The topmost section of data storage is occupied by those variables which appear 
in DIMENSION or EQUIVALENCE statements. The arrangement of this region 
is such that two programs, whose DIMENSION and EQUIVALENCE statements 
are identical, will have this region allocated identically. This fact makes it possible 
to write families of programs which deal with the same data. 

The successively lower sections of storage are occupied by variables not 
mentioned in DIMENSION or EQUIVALENCE statements, then certain storage 
locations required for the operation of the program, and finally a section of 
erasable storage. 

For each object program Fortran produces a printed description of the exact 
arrangement of storage locations. 



Fixed Point 
Arithmetic 



The use of fixed point arithmetic is governed by the following considerations. 

1. Fixed point constants specified in the source program must have magnitudes 
less than 2 15 . 



2. Fixed point data read in by the object program are treated mod 215. 

3. Fixed point arithmetic in the object program is arithmetic mod 2 15 . 

4. Indexing in the object program is mod (size of the object machine). 
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Exponentiation The appearance of an exponential E**F in an expression will produce one of 5 object 

situations. 

1 and 2. If F is a fixed point constant from +1 to +7 inclusive, then one of two 
open subroutines will be compiled, depending upon whether E is fixed or 
floating. One or the other of these subroutines will be compiled for each 
such exponentiation which occurs. These subroutines perform exponentia- 
tion by repeated multiplication; they consist of F— 1 multiplications and a 
maximum of F-f-2 other instructions. The exponential is computed for any 
E. 

3 and 4. If F is a fixed point constant not in the range 4- 1 to -j-7 inclusive, or a fixed 
point variable, then one of two closed subroutines will be compiled, depend- 
ing upon whether E is fixed or floating. These subroutines will not appear 
more than once in the object program. They perform exponentiation by 
forming E F as the appropriate product of E, E 2 , E 4 , E 8 , . . . , followed (if 
F <C 0) by division into I. 

The subroutine for fixed E has 34 instructions and takes an average of 10 
milliseconds. For F ^ 0, E = and for F < 0, lEj > 1, it returns 0; for 
all other cases, it returns E F . 

The subroutine for floating E has 38 instructions and takes an average of 
9 milliseconds. For F <C 0, E = 0, it returns 0; for all other cases, it re- 
turns E F . 

5. If F, and therefore E, is floating, a closed subroutine is compiled. It will 
not appear more than once in the object program. It produces E F as 
£ FinE. j t ft as ioo instructions and takes approximately 6 milliseconds. For 
F <C 0, E = 0, it returns 0; for all other cases it returns |E! F . 



Writing 
Subroutines 
for the 
Master Tape 



Library subroutines exist on the master Fortran tape in relocatable binary form. 
Placing a new subroutine on that tape involves (1) producing the routine in the 
form of relocatable binary cards, and (2) transcribing those cards on to the master 
tape by means of a program furnished for that purpose. 

In the object program transfer to the subroutine is by the sequence 



TSX Subroutine, 4 
Return 



The subroutine itself and any constants that it requires should be located in relocat- 
able locations 0, 1, 2, .... It may also make use of a common storage region of 
any desired length n, beginning with relocatable location 77777 s — (n— 1) and ending 
with relocatable location 77777 8 . 

At the moment of transfer to the subroutine Arg x will have been placed in the 
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Optimisation of 

Arithmetic 

Expressions 



AC, Arg 2 (if it exists) in the MQ, Arg 3 (if it exists) in relocatable location 77775 8 
of the common storage region, Arg 4 (if it exists) in relocatable location 77774 8 , etc. 
The common storage region may also be used for erasable storage by the subroutine. 

The output of the subroutine is to be left in the AC, and index registers 1 and 2 
must be returned with their original contents. 

Fixed point quantities in the object program exist in the following format : sign 
in sign bit, magnitude in decrement field, remainder of word all zeroes. 

It is suggested that error HALTS in subroutines be coded as HPR instructions, 
permitting the tag and address fields to contain identifying numbers which can be 
recognised at the console. 



Considerable attention is given by the Fortran System to the matter of efficiency of 
the object program arising from an arithmetic expression, regardless of how the 
expression has been written. Thus, although, as stated earlier (cf. p. 15 under "Order- 
ing within a Hierarchy"), an expression of the form 



will be taken to mean 



A«B-C'D'E 



■B)-C)-D)-E) 



(where • is * or /, or 4- or — ), it is assumed by the system that any mathematically 
equivalent expression is computationally equivalent. Hence, a sequence of consecu- 
tive multiplications and/or divisions (additions and/or subtractions) not grouped by 
parentheses will automatically be reordered, if necessary, to minimise the number of 
storage accesses in the object program. 

Although the above stated assumption concerning mathematical vs. computa- 
tional equivalence is virtually true in regard to floating point expressions, special care 
must be taken to indicate the order of fixed point multiplication and division, since 
fixed point arithmetic is "greatest integer" (i.e., remain derless) arithmetic in Fortran. 

Thus, the expression 

5*4/2 
is taken to mean ((5*4)/2), but it is computed in a Fortran object program as 

«5/2)*4) 

i.e., from left to right after permutation of the operands to minimise storage accesses. 
The result of the computation in this case is the integer 8. 
On the other hand, the expression 

(5*4)/2 
is taken to mean and, in fact, is computed as 

«5*4)/2) 

and the result is the integer 10. 
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Further Details 
about DO 
Statements 



Thus, to insure maximum accuracy of fixed point multiplication and division with- 
in an otherwise unparenthesised subexpression, judicious insertion of parentheses may 
be indicated. 

Another important type of optimisation, concerned with common subexpres- 
sions, takes place only if the expression has been suitably written. As an example, 
consider the arithmetic formula 

Y = A*B*C+SINF(A*B) 

An efficient object program would form the product A*B only once; yet if the 
arithmetic formula is written as above, the multiplication of A by B will occur twice. 
The correct way to write this arithmetic formula is 

Y = (A*B)*C + SINF(A*B) 

The common subexpression A*B has been displayed by the extra pair of parentheses, 
and an object program will be formed which multiplies A by B only once. 

In general, when common subexpressions exist in an expression, parentheses 
should be used to display them. 

There is one case where the programmer need not write the parentheses, 
because Fortran will understand that they are there. The parentheses discussed in 
the section Hierarchy of Operations in Chapter 3 are of this type, and need not be 

given. Thus 

Y = A*B + C + SINF(A*B) 

is as suitable for optimisation as 

Y = (A*B)+C + SINF(A*B) 

However, the parentheses discussed in the section Ordering within a Hierarchy in 
Chapter 3 must be supplied if common subexpression optimisation is to occur. 

This section contains further details about DOs, which may be of interest to the 
advanced programmer. 

Triangular Indexing. Indexing such as 

DO I = 1,10 
DO J = 1,10 



DO I = 1,10 

DO J = 1,1 
is permitted and simplifies work with triangular arrays. These are simply special 
cases of the fact that an index under control of a DO is available for general use as 
a fixed point variable. 
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Lists for 
Binary Tape 
Operations 



The diagonal elements of an array may be picked out by the following type of 

indexing: 

DO I = 1,10 

Ad, I, I) = some expression 

Status of the Cell Containing I. A DO loop with index I does not affect the contents 
of the object program storage location for I except under certain circumstances, 
namely if 

1. An IF-type or GO TO-type transfer exit occurs from the range of the DO; 

2. I is used as a variable in the range of the DO; or 

3. I is used as a subscript in combination with a relative constant whose value changes 
within the range of the DO. (A relative constant is a subscript the fixed point 
variable of which is not currently under control of a DO.) 

Therefore, if a normal exit occurs from a DO to which cases 2 and 3 do not apply, 
the I cell contains what it did before the DO was encountered. After normal exit 
where 2 or 3 do apply, the I cell contains the first value of the I-sequence which 
exceeds m 2 . After a transfer exit the I cell contains the current value of I. 

What has just been said applies only when I is referred to as a variable. When 
it is referred to as a subscript, I is undefined after any normal exit and is the current 
value after any transfer exit. 

There are restrictions on the complexity of the parenthesised parts of lists for the 
READ TAPE and WRITE TAPE statements which must be observed if the object 
program indexing is not to exceed 288 microseconds, the maximum safe time for 
calculation between CPY instructions on the 704. 

Unfortunately this matter is exceedingly complicated and therefore no complete 
discussion will be given. Instead, certain rules will be given which will permit the 
construction of lists which can be relied upon to work; there will be other lists, 
excluded by these rules, which would also work. 

Define the term subscript combination as follows. Consider each subscript of 

each subscripted variable to be specified by the four symbols, c, v, c' and d, where 

d is the maximum value of the subscript as given in the DIMENSION statement 

about the variable, and where c*v±c' is the full form of the subscript. (The full 

form of the subscript I would be 1*1 + 0; that of the subscript 3 would be 1*0 + 3, 

etc.). Then the subscript combination of a subscripted variable is one of the 

ordered sets 

Ci, v, 

Ci, v 1( c 2 , v 2 , d! 

Ci, v b c 2 , v 2 , c 3 , v 3 , d,, d 2 

depending upon whether the variable has 1, 2, or 3 subscripts. Thus A(I,J,K), 
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B(I,K,J), C(5*I,J,K), and D(I,J) all have different subscript combinations, but 
A(I,J,K) and B(I+5,J,K) have the same subscript combination if the first two 
dimensions of A and B are the same. 

Also, define as an element of a list a part bounded by principal commas in 
the list. 

Elements which Contain no Subscripted Variable. There are no restrictions on 
such elements. 

Elements which Contain no Variable with Three Subscripts. Let 

N, = number of different 1-dimensional subscript combinations 

N 2 = " " " 2- 

Mt = if N, = 0; otherwise M, = 1 

M 2 = if N 2 = 0; otherwise M 2 = 5 

Then the restriction is that 4N- -!-6Nf -i-max^M- M^ must not exceed 18. 

For example, the element ((A(I,J), B(I,J), J = 1,10), C(I), I = 1,10) has 
N x = 1, N 2 = 1 (provided that the first dimensions of A and B are the same); 
hence 4N 1 + 6N 2 +max(M 1 ,M 2 ) = 15, and the element is permissible. 

Elements which Contain Variables with Three Subscripts. Consider a variable A 
whose three subscripts s l9 s 2 , s 3 involve as variables v l5 v 2 , v 3 . Then the skeleton 

elements 

( ( ( , v, = ), v 2 = ), , v, = ) 

or 

( ( ( , V 3 = ), V, = ), , v 2 = ) 

or 

( ( ( , V 2 = ), V 3 = ), , V, = ) 

will work. The innermost parentheses (controlling v 3 , v 3 , v 2 respectively) may 
contain A(s 1 ,s 2 ,s 3 ) and any number of other variables with the same subscript 
combination; similarly the outermost parentheses may contain variables of any one 
subscript combination. The middle parentheses may not contain any subscripted 
variables. Thus for example 

( ( (A(I,J,K), B(I,J,K), K = 1,10), J = 1,5), CO), D(l), I = 1,10) 

will work provided that the first two dimensions of A and B are the same. This 
example makes use of the first of the skeleton elements just given. 

In the special case where the next element in the list is neither subscripted nor 
enclosed in controlling parentheses, the first of the above skeleton elements may 
also contain variables of any one subscript combination in its middle parentheses. 

Finally, the skeleton element 

( ( ( , v, = 1,(1,), v 2 = l,d 2 ), v 3 = ) 
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which unlike the others indexes in the natural order, will work. Notice, however, 
that the indexing parameters for v 1 and v 2 must be such that the array is swept 
through consecutively. The innermost parameters may contain A(s l5 s 2 ,s 3 ) and any 
number of other variables with the same subscript combination. 

Variable Indexing Parameters and Relative Constants. Another restriction affect- 
ing lists for binary tape operations concerns the use of relative constants and of 
indexing parameters which are variables. (A relative constant is a subscript, the 
variable of which is not currently under the control- of a DO or a controlling 
parenthesis. ) 

The restriction is, that variables which have a subscript involving either a 
relative constant or an index governed by variable indexing parameters may appear 
only in the first element of a list. Furthermore, all such variables must also have 
among their subscripts one whose index is controlled by the first parenthesis. 

For example, in ((A(I,J), I = 1,L), B(K,J), J = 1,M) both the subscripts of 
A are governed by variable indexing parameters; for B the same is true of one of its 
subscripts, while the other is a relative constant (unless the READ TAPE or WRITE 
TAPE is itself in the range of a DO for K). However, since J, which is the index 
controlled by the first parenthesis, is a subscript of both A and B, this element will 
work if it is the first element in the list. 



Concept of 
Relative 
Constant 
Definition 



If a fixed-point variable appears in a relative constant anywhere in the program, then 
any appearance it makes as the left-hand side of an arithmetic formula or in an input 
list will be called a relative constant definition. In the object program some computa- 
tion will take place at each such definition, and if this definition occurs in a READ 
TAPE or READ DRUM list, a timing violation may occur unless the following 
restrictions are observed. 



Restrictions 
on Relative 
Constant 
Definitions 
in READ TAPE 
Lists 



There are two general cases to be distinguished. 

Case 1. The relative constant occurs in the program only elsewhere than later in 

the same list with the relative constant definition. 
Case 2 . The relative constant occurs in the same list with the relative constant 

definition and following it. 

Case 1. Relative Constant Elsewhere than Later in the List. There are two subcases. 
1. Where the relative constant occurs in a 2- or 3-dimensional subscript combination, 
or in a 1 -dimensional subscript expression with a coefficient, the following rule 
must be observed. 

Rule. The relative constant definition may be in a READ TAPE list, in any 
position, except (a) to the left of any DO-implying parenthesis (i.e., a left paren- 
thesis which is not a subscript parenthesis — see page 26), (b) to the left of a 
non-subscripted array symbol, or (c) within DO-implying parentheses. 



47 



2. Where the relative constant occurs in a 1 -dimensional subscript expression without 
a coefficient, the following rule must be observed. 

Rule. Relative constant definitions in a READ TAPE list are permissible, but 
there must not be more than five of them where the conditions of the three 
exceptions above apply. 

Case 2. Relative Constant Later in the List. The relative constant may follow its 
definition in the same list only under the following special conditions. 

1. The relative constant appears only later in the same list and not elsewhere in the 
program. 

2. The relative constant is 1-dimensional and without a coefficient. 

3. An extra pair of parentheses is placed around the relative constant and the symbol 
which it subscripts. 

Examples: The READ TAPE list I, J, K, (B(J, L)) is illegal because the relative 
constant is 2-dimensional. The list I, A (I), is illegal but becomes legal with extra 
parentheses; i.e. I, (A (I) ). 



Restriction 

on Relative 

Constant 

Definitions 

in READ DRUM 

Lists 



Relative constant definitions may not be to the left of an array symbol in a READ 
DRUM list. 



Restrictions 
on Relative 
Constant 
Definitions 
in READ and 
READ INPUT 
TAPE Lists 



Where the relative constant definition appears in the same READ or READ INPUT 
TAPE list with its relative constant and precedes it, extra parentheses may be re- 
quired in the list. In such a list, it is necessary that there be a left parenthesis, other 
than the left parenthesis of a subscript combination, between the relative constant 
definition and its relative constant. If the list does not contain the parenthesis, it may 
most easily be obtained by placing parentheses around the symbol subscripted by 
the relative constant. 
Examples : 

A, B, K, M, (C(J), J = 1, 10), G(K) 
A, B, K, M, GOO 

The first of these two BCD input lists is correct. The second is incorrect, but may be 
made correct with extra parentheses; i.e. 

A, B, K, M, (G(K)) 



Limits on the 
Size of the 
Source Program 



In performing the translation from source to object program, Fortran forms and 
uses tables which summarise various aspects of the information contained in the 
source program. These tables are limited in size, with corresponding limitations on 
the amount of information which the source program may contain. If a table size is 
exceeded, the Fortran program will HALT at a stop location. 
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In what follows, the phrase "literal appearance" means that if the same thing appears more 
than once it must be counted more than once. 

1. (TEIFNO Table). The number of Fortran statements which have statement numbers 
must not exced 1500. (An input or output statement which has a statement number 
and whose list contains controlling parentheses counts double.) 

2. (FIXCON Table). The number of different fixed point constants must not exceed 100. 
(In this count, constants which differ only in sign are not considered different.) 

3. (FLOCON Table). The number of different floating point constants must not exceed 
450. (In this count, constants which differ only in sign are not considered different, nor 
are numbers such as 4., 4.0, 40.E-1, etc., which are really the same number.) 

4. (TDO Table). The total number of DOs in a program must not exceed 150. It must be 
remembered that a DO derives not only from a DO statement but also from DO-im- 
plying parentheses in input-output lists. 

5. The number of DO statements in any one nest must not exceed 50. 

6. (TIFGO Table). The total number of ASSIGNS plus IF-type and GO TO-type state- 
ments must not exceed 300. 

7. (TRAD Table). The total number of statement numbers mentioned in assigned GO TO 
and computed GO TO statements must not exceed 250. 

8. (FRET Table). The total number of numbers mentioned in FREQUENCY statements 
must not exceed 750. (Such a statement as FREQUENCY 30(1, 2, 1) has 4 numbers.) 

9. (DIM Tables). The total number of 1-, 2-, and 3-dimensional variables which appear 
in DIMENSION statements must not exceed 100, 100, and 90, respectively. 

10. (EQUIT Table). The total number of literal appearances of variables in EQUIVA- 
LENCE statements must not exceed 750. 

11. (FORMAT Table). Compute the value f in the following way for each FORMAT 
statement in the problem. Make a character count of each FORMAT statement, start- 
ing at the opening parenthesis and ending at the closing parenthesis. Divide by 6, 
allowing a remainder to count as 1 . Add 1 . The total f values for each problem must 
not exceed 715. 

12. (LAMBDA Table). This table, and the BETA table which follows, limit the size of 
arithmetic expressions on the right-hand side of arithmetic formulas and as the argu- 
ments of IF statements. In any one expression, let 

n = number of literal appearances of variables and constants, except those in subscripts; 

b = number of open parentheses, except those introducing subscripts; 

p = number of appearances of + or — , except in subscripts or as unitary operators. 

y A 111^ "J - JX1 1 \ "'■ V. [ *-* ) U U lAillLUi J* Oj-'WltiLV^J. - J , 

t = number of appearances of * or /, except in subscripts; 

e = number of appearances of * * ; 

f = number of literal appearances of function names; 

a = number of arguments of functions. (For SINF(SINF(X)), a = 2.) 

Then X, which equals n+4b-f-4a — 3f+3p-f 2t+e+3, must not exceed 400. 

13. (BETA Table). With the same definition, /? = X+l— n — f must not exceed 300. 
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14. (ALPHA Table). In order to determine in advance whether an ALPHA table overflow 
will occur during the course of translation of an arithmetic statement, the following 
procedure should be carried out. Set the initial value of a counter to 3. Scanning the 
right-hand side of the statement in question, add 4 to the value of this counter for each 
left parenthesis encountered and subtract 4 for each right parenthesis encountered. This 
statement is compilable by the 4K version of the Fortran System if and only if the 
counter value never exceeds 43. This statement is compilable by the 8K version of the 
Fortran System if and only if the counter value never exceeds 139. 

15. (CLOSUB Table). In the entire program, the number of literal appearances of func- 
tions must not exceed 1500. 

i6. (FORVAL Table). The total number of literal appearances of non-subscripted fixed 
point variables on the left-hand side of arithmetic formulas and in input lists must not 
exceed 500. 

17. (FORVAR Table). The total number of literal appearances of non-subscripted fixed 
point variables on the right-hand side of arithmetic formulas and in the arguments of 
IFs must not exceed 750. 

IS. (FORTAG Table). The iotai number of literal appearances of subscripted variables 
must not exceed 1500. 

19. (TAU Tables). The total number of different 1-, 2-, and 3-dimensional subscript 
combinations must not exceed 100, 90, and 75, respectively. (See the section "Lists for 
Binary Tape Operations" for the definition of subscript combination.) 

20. (SIGMA Tables). Consider a variable with 3 subscripts and let the additive "arts 
of these subscripts (when written in full form) be c/, c 2 ', c 3 '. Then the number of 
distinct ordered triples (c/, c/, c/), such that at least one c/ does not equal zero, must 
not exceed 99. Similarly, for 2- and 1 -dimensional variables, neither the number of 
distinct ordered couples (c/, c 2 ') nor of distinct c/ may exceed 99. 

21. (FORSUB Table). The total number of distinct Fortran functions (functions defined 
within the problem by an arithmetic statement) must not exceed 50. 

22. (NLIST Table). The total number of unique fixed point variables in assigned GO TO 
statements must not exceed 25. 

The following are tables applying exclusively to Fortran II: 

23. (SUBDEF Table). This table arises from arguments given in SUBROUTINE and 
FUNCTION statements. One word is used for each SUBROUTINE or FUNCTION 
name and one word for each literal appearance of an argument. The number of words 
must not exceed 900. 

24. (COMMON Table). The number of literal appearances of variables in COMMON 
statements must not exceed 600. 

25. (HOLARG Table). This table arises from Hollerith arguments in CALL statements. 
For every nH in a CALL statement, divide n by 6, adding 1 for any remainder. Add 
1 to this. This total for any one program must not exceed 900. 

26. (NONEXC Table). The total number of non-executable statements must not be greater 
than 1500. 

27. (TSTOPS Table). The total number of STOP statements must not be greater than 
300. 
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CHAPTER 8. EXAMPLES OF FORTRAN PROGRAMMING 



A Complete 
but Simple 
Program 



The example below is the same as that given in Chapter 1 , but expanded into a 
complete program. The purpose of the program is to discover the largest value 
attained by a set of numbers A (I) and to print that number on the attached 
printer. The numbers A (I) exist on punched cards, 12 to a card, each number 
occupying a field of 6 columns. There are not more than 999 numbers; the 
actual number N is punched on a lead card. 



a FOR 

c-* 

COMMENT 


| 


FORTRAN STATEMENT 


IDENTI- 


STATEMENT 


7 


73 HI 


C| 




PROGRAM FOR FINDING THE LARGEST VALUE 




cj 


X 


ATTAINED BY A SET OF NUMBERS 








DIMENSION A(999) 




I 




FREQUENCY 30(2,1,10), 5(100) 








READ 1, N, (A(I), I - 1,N) 




J 1_ 




FORMAT (I3/U2F6.2)) 




I BIGA - A(l) 




! 5 


! DO 20 I = 2 , N 




| 30 


1 IF (BIGA-A(I)) 10.20,20 




1 10 '■ BIGA - Ad) 




1 

| 20 




CONTINUE 




1 


PRINT 2, N, BIGA i 


1 2 


FORMAT (22H1THE LARGEST OF THESE 13, 12H NUMBERS IS F7.2? 


J I STOP 77777 




1 : i „.,._ . 





The first executable statement is the READ; therefore the program begins there. 

The READ causes first N and then A(l), A(2), , A(N) to be brought 

in from the card reader, in accordance with the FORMAT statement 1 . Notice 
that the indexing of the loop bringing in A will work corr^*y, since the 
indexing parameter N occurs earlier in the list. 

The FORMAT statement says that there is first a single card with format 
13, followed by any number of cards with format 12F6.2. On the single card 
N is punched as a decimal integer in columns 1-3 and is also to be converted 
into a fixed point number. In the remainder of the deck the numbers A (I) are 
punched 12 to a card in columns 1-6, 7-12, etc. Each number is presumably 
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punched as xxxxx or — xxxxx with the decimal point understood to precede the 
last two digits. If, however, a column is used for a decimal point, its position 
overrides this understanding. 

After the READ is executed control moves to the next executable state- 
ment, the arithmetic formula BIGA = A(l). The cell BIGA now contains A(l). 

Next the DO statement sets I to 2 and creates a loop starting with the IF 
and ending with the CONTINUE. The first time the IF is executed it transfers 
control either to statement 10 or statement 20, according as A(2) > BIGA or 
A(2)^ BIGA. In the first case, therefore, BIGA becomes A(2); otherwise it 
remains A(l). 

Control is now in the CONTINUE, which is the last statement in the range 
of a DO which is still unsatisfied (provided that N > 2), and which is not a 
transfer. Therefore DO-sequencing occurs : I is increased to 3 and control goes 
back to the IF. The comparison now is between BIGA and A(3); ifA(3) is 
the larger it becomes the new BIGA. 

This process continues until control is in the CONTINUE with I = N. 
The DO is satisfied; therefore normal sequencing occurs and control moves to 
the PRINT. BIGA now contains the largest value in the set of A's; N and the 
A-array have not been altered. 

The PRINT statement causes N and BIGA to be printed on the attached 
printer in accordance with the FORMAT statement 2. This statement causes 
a line to be printed as follows. The carriage control character is a Hollerith "1" 
which causes a skip to channel 1, bringing the paper to the top of a new sheet. 
Type wheels 1-21 receive the text 'THE LARGEST OF THESE '; wheels 22-24 
receive N converted back to decimal integer form; wheels 25-36 receive ' NUM- 
BERS IS '; and wheels 37-43 receive BIGA converted into fixed decimal form, 
with a decimal point preceding the last two digits. Notice that an extra column 
is allotted to BIGA to allow for the decimal point, which was presumably 
omitted in the data input. 

The program ends with a HALT and 11111$ in the address field of the 
storage register on the console. Pressing the START button will have no effect. 

The only subscripted variable used in the program is A. Hence DIMEN- 
SION A (999) is the only dimension information required. 

The only statements about which a FREQUENCY statement can be made 
are the IF and the DO. The programmer anticipates that BIGA will usually be 
greater than A (I), a plausible guess if the A (I) are randomly arranged but 
not if they tend to increase with increasing I. He also predicts for the DO that 
100 is a reasonable average for N. 
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A DO Nest 
with Exit 
and Return 



Given an N x N square matrix A, to find those off-diagonal elements which are 
symmetric and to write them on binary tape. 



_ FOR 

COMMENT 


o 

o 
6 


FORTRAN STATEMENT 

7 72 


IDENTI- 
FICATION 

73 SO 


STATEMENT 
NUMBER 

1 5 


1 




REWIND 3 








DO 3 I = 1 , N 






— 


DO 3 J ~ 1 , N 






IF(A(I,J)-A(J, I)) 3,20,3 




1 3 




CONTINUE 








END FILE 3 








MORE PROGRAM 












| 20 




IF(I-J) 21,3,21 




I 21 




WRITE TAPE 3, I, J. A(I,J) 






— 


GO TO 3 




...i. 







After rewinding tape 3, a nested pair of DO loops scans the entire matrix 
for elements A(I,J) equal to A(J,I). Whenever such an element is found an 
exit completely out of the nest is made to a routine which for off-diagonal 
elements only writes a 3-word record (I, J, and A(I,J)) in binary on tape 3. 
Both for on- and off-diagonal elements this routine makes no change in the 
indexes or indexing parameters of the nest, and so it is permissible to re-enter 
the nest and continue the scan. 

This program actually finds each element twice. This could be avoided by 
writing the second DO as DO 3 J = I,N. 
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APPENDIX A. TABLE OF FORTRAN CHARACTERS 





0. 






Q. 






a 






Q. 






4 






4 






4 






4 







H 




a 


H 




a 


»- 




a 


H 




tt 


a 


« 


B 


a 


* 


K 


a 


<t 


K 


a 


<t 


4 


o 


o 


4 


o 


o 


4 


o 


o 


< 


o 


e 


o 


a 


i* 





n 


r» 


O 


m 


r« 


o 


ffl 


N 



1 


1 


01 


01 


A 


12 

1 


61 


21 


J 


11 
1 


41 


41 


/ 




1 


21 


61 


2 


2 


02 


02 


B 


12 
2 


62 


22 


K 


11 
2 


42 


42 


S 



2 


22 


62 


3 


3 


03 


03 


C 


12 
3 


63 


23 


L 


11 
3 


43 


43 


T 



3 


23 


63 


4 


4 


04 


04 


D 


12 

4 


64 


24 


M 


11 
4 


44 


44 


U 




4 


24 


x 64 


5 


5 


05 


05 


E 


12 
5 


65 


25 


N 


11 
5 


45 


45 


V 



5 


25 


65 


6 


6 


06 


06 


F 


12 
6 


66 


26 





11 
6 


46 


46 


W 



6 


26 


66 


7 


7 


07 


07 


G 


12 
7 


67 


27 


P 


11 
7 


47 


47 


X 



7 


27 


67 


8 


8 


10 


10 


H 


12 
8 


70 


30 


Q 


11 
8 


50 


50 


Y 



8 


30 


70 


9 


9 


11 


11 


1 


12 
9 


71 


31 


R 


11 
9 


51 


51 


Z 



9 


31 


71 


Blank 




20 


60 


+ 


12 


60 


20 


- 


11 


40 


40 








12 


00 


= 


8-3 


13 


13 




12 
8-3 


73 


33 


$ 


11 
8-3 


53 


53 


- 



8-3 


33 


73 


- 


8-4 


14 


14 


) 


12 

8-4 


74 


34 


* 


11 
8-4 


54 


54 


( 



8-4 


34 


74 



Note 1. There are two — signs. Only the 11 -punch may be used in source 
program cards. Either minus may be used in input data to the object program; 
object program output has the 1 1 minus. 

Note 2. The $ character can be used in Fortran only as Hollerith text in a 
FORMAT statement. 
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APPENDIX B. TABLE OF FORTRAN STATEMENTS 



STATEMENT 


NORMAL SEQUENCING 


a = b 


Next executable statement 


GO TO n 


Statement n 


GO TO n, (n,,n 2 ,...,n m ) 


Statement last assigned 


ASSIGN i TO n 


Next executable statement 


GO T0(n lr n 2 ,...,n m ), i 


Statement nj 


IF (a) n 1r n 2 ,n 3 


Statement ni,n 2 ,n 3 as a less than, =, or greater than 


SENSE LIGHT i 


Next executable statement 


IF (SENSE LIGHT i) n„n 2 


Statement n 1; n 2 as Sense Light i ON or OFF 


IF (SENSE SWITCH i) n,,n 2 


" " as Sense Switch i DOWN or UP 


IF ACCUMULATOR OVERFLOW n,,n 2 


" " " as Accumulator Overflow trigger ON or OFF 


IF QUOTIENT OVERFLOW n,,n 2 


" " as MQ Overflow trigger ON or OFF 


IF DIVIDE CHECK n,,n 2 


" " as Divide Check trigger ON or OFF 


PAUSE or PAUSE n 


Next executable statement 


STOP or STOP n 


Terminates program 


DO n i = m 1 ,m 2 or DO n i = mi,m 2 ,m 3 


Next executable statement 


CONTINUE 


FORMAT (Specification) 


Not executed 


READ n, List 


Next executable statement 


READ INPUT TAPE i, n, List 




PUNCH n, List 




PRINT n, List 




WRITE OUTPUT TAPE i, n, List 




READ TAPE i, List 




READ DRUM i, j, List 




WRITE TAPE i, List 




WRITE DRUM i, j, List 




END FILE i 




REWIND i " 




BACKSPACE i 




DIMENSION v,v,v 


Not executed 


EQUIVALENCE (a,b,c,.-), (d.e.f,..), ... 


FREQUENCY n(i,j,.. ), m(k,l,..), ... 
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INDEX 



Admissible characters 54 
Arithmetic iormuias il, iu, 38 

maximum length of 49 

mode of result 16 

not equivalent to ASSIGN 18 

optimisation of 43 
Arrays 9, 10 

arrangement in storage 1 1 

order of reading or writing 

25, 26, 34, 46, 47 

size of 11,37 
Assigned GO TO in range of DO 

18. 23 

Blank characters 7, 8, 29, 30 

Card punching 7 
Carriage control 30, 49 
Checking of reading 31, 32, 33 

of writing 32 
Comment cards 8 
Common storage region 42 
Common subexpressions 44 
Constants 9 

maximum number of 49 
Continuation cards 7 

Data input to object program 

30,51,54 
DO-sequencing 24 

examples of 7,21,23,51,53 
DO statements 21,44 

and transfers 22, 24 

exit and return 22, 53 

index of 21,23,24,26, 

38, 39, 40, 44 

maximum number of 49 

nest of 22, 49, 53 

normal exit from 23, 45 

permitted calculation in 

range of 23 

preservation of index 23,39,45 

range of 21, 22 



satisfied 21, 24, 52 
transfer exit from 22, 23, 45 
triangular indexing 44, 53 
within DOs 22 

Echo checking 32 
End of File 8,31,32,33,35 
Exponentiation 9, 15, 28, 30, 42 
Expressions 14. 19 

maximum length of 49 

mixed 14, 16 

optimisation of 43 

rules for forming 14 

Format 27 

exam n Ies of 53 52 

multi-record 29 

of input data to object program 

30,51,54 

of fixed point quantities in 

storage 43 
Function statements 12, 17 
Functions 3, 12 

built-in 13 

maximum number of 50 

modes of 12 

rules for writing 42 

HALTS, identification of 

20,21,36,43,48,50 
Hierarchy of operations 

15, 43, 44 
Hollerith fields 29, 52, 54 

Indexing parameters 22, 23 

26, 39,40,48,51,52 
Input to Fortran 8, 54 

Library subroutines (see 
Functions) 

Lists of quantities for 
transmission 25, 33, 34 
abbreviated form 26, 34 



for binary tape 27, 33, 34, 

3C AC 
J, tj 

for drum 27, 34, 35 

Non-executed statements 

23, 24, 27, 37, 50, 55 
Normal sequencing 24, 55 

Object program arrangement 
in storage 41 
Order of arrays 1 1 

of hierarchy of operations 15 

of operations within 

hierarchy 15 

of quantities in lists 25, 26, 

34, 45, 46, 47 

of statements 8 

Parentheses 14, 15, 26, 28, 29, 

44, 47, 48 
Precision 3, 9, 10, 41 

Redundancy checking 32, 34 
Relative constants 45, 47, 48 

Scale factors in input-output 28 
Sequencing, summary of 24, 55 
Size of machine required 2,41 
"Small" 704 2, 41 
Statement numbers 8, 9 
Statements 7, 8 

maximum number of 49 

summary of 8, 55 
Storage allocation 11, 37, 38, 41 
Subscript combination 45,46,47 

maximum number of 50 
Subscripts 9, 10, 14, 46, 50 

maximum value of 11,37,50 

minimum value of 1 1 

Table sizes 47, 48 

Variables 10 

maximum number of 49, 50 
restrictions on naming 10 
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